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CNI | In this paper we show that if the refractive index, or rather, [n(oj) — 1] satisfies the dispersion relations 
then, it is implied by Titchmarsh's theorem that, n(uj) — > 1 as u) — > oo. Any other limiting value for n(io) 
would violate relativistic causality, by which we mean not only that cause must precede effect but also that 
signals cannot travel faster-than-c, the velocity of light in a vacuum. 
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q ■ 1 Introduction 

o: 

r—i There has been much interest in recent theory and experiment on superluminal light propagation in media 
CO [1, 2, 3, 4, 5, 6, 7]. These articles appear not to violate relativistic causality, by which we mean the front 
. velocity cannot travel faster-than-c. The front velocity is the speed at which the very first extremely 
small vibrations of the wave will occur [8] . The front of the wave is composed of the highest uv frequency 
Oh, components (ie u> — + oo). Sommerfeld [9], predicted that the front velocity should travel at c. A brief history 
of events is given in Brillouin's book [8]. Later Voigt [8], gave a physical explanation for Sommcrfelds result, 
stating that the modern theory of dispersion uses the assumption of point like electrons with a finite mass. 
' Inertia prevents the electrons from oscillating at the very start of the wave. Electron oscillation needs time 
£7^. to build; only after the wave has been in motion for some time can the electrons react back on the wave. The 
very early oscillations of the wave therefore pass though the dispersive medium as if through vacuum. Also, 
electrons having a finite mass, cannot oscillate at infinite frequencies without amassing infinite energy in the 
process. We will restrict our material media to less than infinite energy. If n(uo) — > (3 as ui — > oo we take the 
front velocity to be v — c/{3, it is required that (3=1. In fact we will prove (via Titchmarsh's theorem to be 
stated later) that in order that n(u>) — 1 not violate causality it must obey the dispersion (Kramers Kronig) 
relations and as a consequence of this n(uj) — > 1 as u> — > oo. 

This may appear as obvious and well known to some readers. However several well known text book 
accounts [10, 11, 12] on the subject invoke a physical model for refractive index to require n(u>) — ► 1 at 
high frequencies as we shall discuss later. It has not been made clear that the mathematical derivation of 
the dispersion relation for [n(u>) — 1] requires [n(oj) — 1] — > as lo — > oo as part of the proof. The precise 
physical model is irrelevant. This is an important point which this article will clarify. 

It appears not all recent articles agree that n(o>) — > 1 for high frequency is a requirement for causality of 
n(io) — 1. It has been claimed that a signal can travel at v = c/n(oo) where n(oo) < 1, [13]. Furthermore, 
it was first suggested that the Scharnhorst effect [14, 15, 16, 17], relating to the normal propagation of light 
between two parallel mirrors in a vacuum might allow a signal velocity (more precisely a front velocity) to 
exceed c. In the upcoming section we briefly discuss the Scharnhorst effect and give arguments which show 
that the original derivation may be flawed. Also we show that the effect is extremely small and in practice 
unmcasurable. We then concentrate on physical media (non-vacuum) where experimental verification is 
possible. Our purpose is to show that in general signals travel at c for material media despite experimental 
claims to the contrary in the literature. We give physical arguments to show why u(lo) — > 1 as lo — ► oo for 
these physical media. 
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2 Scharnhorst Effect 



The Scharnhorst effect relates to light propagating in the vacuum between two parallel mirrors. The vacuum 
modes are changed by the boundaries (in much the same way as in the Casimir effect) and the light expe- 
riences the vacuum as a dispersive bircfringcnt medium. The real part of the refractive index in a direction 
parallel to the mirror surface is unity. The refractive index in a direction perpendicular to the mirror surface 
is found to be less than unity. These calculations are done using perturbation theory valid only for small 
frequencies uj « m where m is the electron mass. 

Scharnhorst [14], derived a refractive index (perpendicular to the mirror surface) for the vacuum based 
on low frequencies and showed n(0) < 1. (Note that n(0) implies the small frequency limit.) Combining 
this with the Kramers Kronig relations written for Re n(ui)— Re n(oo) (as opposed to Re [n(u>) — 1]), and 
setting u> = 0, Barton and Scharnhorst [16] showed that either (real part) Re n(oo) < 1, which would imply 
signals moving at faster-than-c speed, or (imaginary part) Im n(u>) < for some frequency range, which 
would imply that the vacuum could amplify the light signal for some range of frequency. The paper [16], 
discusses both options but does not make a definite choice between the two. 

The dispersion relations used by Barton and Scharnhorst [16] and later by Scharnhorst alone [17] are not 
complete in the sense that they have set a term Im n(oo) = 0, which in principle might change their end 
result. (See Eq. (41) in this paper). Barton and Scharnhorst say [16]; 

"We have no conclusive proof that Re n converges and Im n vanishes as ui — > oo because [• • •] 
their asymptotics are quite likely to be governed by non-perturbative effects that we cannot 
calculate." 

The reasons for assuming Im n(w — > oo) = are given in their paper [16] and we refer the reader there 
for further discussion. A further objection to the Scharnhorst effect is that it violates special relativity (SR). 
Even if you cannot measure the effect, to be discussed below, the fact that it can exist even in principle is 
objectionable. 

Consider the case of a light clock, a pair of mirrors with a light pulse bouncing between them. This type 
of clock has been used to derive effects like time dilation and Lorentz contraction in undergraduate text 
book accounts of SR. For time dilation the clock moves in a direction parallel to the surface of the mirrors. 
Since the Scharnhorst effect predicts that the vacuum refractive index is equal to unity in that direction 
then we would expect the light clock (in a Scharnhorst type arrangement) to give much the same prediction 
as in SR. It should be noted however that the speed of light in the moving frame of the clock has altered 
due to the Scharnhorst effect, so predictions are not exactly the same as with no Scharnhorst effect. In the 
Lorentz contraction arrangement, the light clock is tilted on its side. The direction of motion of the clock 
is the same as the light bouncing between the mirrors, (ie. perpendicular to the mirror surface). It appears 
that several light clocks in relative motion would not calculate the same length scales if the Scharnhorst 
effect were operating between the mirrors. Observers in relative motion would not know what to use for c 
in the SR velocity addition formula. The change in the velocity of light, as predicted by the Scharnhorst 
effect, is inversely proportional to the fourth power of the distance between the mirrors, see 5c below. The 
distance measured between the mirrors differs from one reference frame to another in relative motion due to 
Lorentz contraction. Hence, the measured value of the velocity of light, measured from one frame to another, 
must also change and this violates SR in a very fundamental way. The whole of SR is based on c being an 
invariant. Which invariant value of c do you use in this case? It appears SR does indeed break down if the 
Scharnhorst effect were real, even if you cannot measure the Scharnhorst change in c in practice. 

Two papers have appeared both stating that measurement of faster-than-c signals between mirrors was 
impossible. The first was by Milonni and Svozil [18], which uses an argument based on the uncertainty 
relation for velocity and the uncertainty in time due to switching on a signal. Their argument can be 
summarized as follows; A measurement of velocity involves a distance and time measurement v = L/t. The 
time is limited by a signal turn on time St w l/u> where uj is the frequency of the signal. We assume the 
shortest possible delay. Then Sv = L5t/t 2 > c 2 5t/L. Using c/co ~ A we obtain Sv = cX/L. The change in the 
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velocity of light predicted by the Scharnhorst effect is 8c — kca 2 (A c / L) 4 where \ c = %/ (rac) = 3.9 x 10 11 cm 
is the Compton wavelength, a = 1/137 and k w 10~ 2 . Hence the ratio of 

TT>T^h- =1.5xlO« - (1) 



<5c /ca 2 \A C / \A 

where we have used A w A c . Thus the measured uncertainty in velocity, <5w, is much greater than the 
predicted change in the signal velocity of light, 5c, by many orders of magnitude. Thus the Scharnhorst 
effect is not useful as an experimental verification that faster-than-c speeds for signals are possible. Milonni 
and Svozil also conclude [18], 

"... it is clear that the uncertainty in the measured propagation velocity will always be enor- 
mously larger than the correction to c associated with the Scharnhorst effect. We conclude, 
therefore, that no measurement of the faster-than-c velocity of light predicted by the Scharn- 
horst effect is possible." 

The second paper by Ben-Menahem [19], uses an argument based on the sharpness of the wavefront. 

"... in order to observe faster-than-c propagation of the wavefront, it is necessary to sharpen 
the falloff of the fields at the wavefront to a length scale less than 1/m [where m is the electron 
mass]. This feat requires the inclusion in the packet of waves with u> > m, for which eq. (3) [the 
vacuum refractive index derived by Scharnhorst for uo <C m] is a bad approximation." 

In the following sections, we will concentrate on physical media (non-vacuum media which have only finite 
energies), in flat space-time, and we will derive the dispersion relations for n(u>) — 1 and in the process prove 
that for causal propagation n(co) — > 1 asw^oo. We will state Titchmarsh's theorem and use it to prove 
our result. We conclude with a discussion that for non-violation of causality it is a necessary condition that 
for physical media n(co) — > 1 for lo — > oo and thus signals must travel at c, the velocity of light in free space. 



3 Pulse propagation through a medium 

Consider an electric pulse E(z,t) incident on a thin dielectric medium of thickness 5. The following argu- 
ments can be found in the book by Nussenzveig [20]. The pulse is normally incident on the medium and 
has a transmission coefficient given by 2/(n(u>) + 1) which we shall take to be approximately equal to unity. 
This turns out to have no influence on causality considerations. We assume that the medium has a complex 
refractive index n(u) = n r {u) + irn(uj). We define the absorption coefficient a = 2rii(uj)uj / c. 

1 f°° 

E(z,t) = —j E in (z,u;)e- Mt cLo (2) 

The pulse is made up of many Fourier components. The component of frequency to can be written as 

£ in (z,^H£ in (0,^) e ^/ c (3) 

If the pulse travels along the z axis and we place the left face of the medium at position z = then the 
corresponding electric field component on the exit face of the medium is 

E out (S,u J ) = E iQ (0,Lo)e^ n ^ S / c (4) 

We will assume that the output is a linear function of the input, so that the superposition principle applies. 
This shows that given input f(t) the output x(t) must satisfy 



/oo 
g(t,t')f(t')dt' (5) 
-oo 

Furthermore we assume time translation invariance so that if the input is delayed (or advanced) by a time 
interval t then the output will be delayed (or advanced) by the same time interval. It follows that the function 
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g(t, t') must depend only on the time difference between the input and output pulses, g(t, t') = g(t — t'), so 
that 

/>OC 

x(t) = / g(t - t')f(t')dt> (6) 



This takes the form of a familiar convolution integral ("Faltung" integral in German, meaning folding inte- 
gral). Equation (6) has the consequence that 

X{w) = G(u)F(uj) (7) 

where F(u)),X(u) and G(w) are the Fourier transforms of f(t),x(t) and g(t) respectively. Thus 

/oo 
g(r)e iUJr dr (8) 
-oo 

If we invoke a primitive causality condition; that the output cannot precede the input, so that if f(t) = 
for t < t then the same is true for x(t). Hence g{r) — for r < and we can replace the lower bound in 
the above integral by zero. 

/>oo 

G(uj) = / g(r)e^ T dr (9) 
Jo 

This implies that G(u>) has a regular analytic continuation over the upper half of the complex plane, i+. If 
we write w = x + iy for y > then 

/.oo 

G(w) = / g{r)e ixr e~ VT dr (10) 
Jo 

the decay factor exp(— yr) can only help convergence when r > 0. We still do not have sufficient conditions 
to derive the dispersion relations for G(ui). We need information regarding the limiting form of G(uj) as 
u> — > oo. In particular, we need to prove that it decreases rapidly at infinity, otherwise the Re G(u>) and Im 
G(to) may be totally unrelated as in the case of a complex constant a+ib, [20]. We take G(co) to be square 
integrable along the real axis, (in mathematical notation, G(co) G L 2 (— oo,oo)), so that 



f 

J —i 



\G{lo)\ 2 (Lj < K (11) 



where if is a finite constant. If F(uj) and G(u) are Fourier Transforms of f(t) and g(t) defined in a similar 
fashion to E- m (t) then according to Parseval's theorem 

/OO -| /.oo 

f(t)g*(t)dt = — / F(w)G*(uj)dw . (12) 
-oo ^ J —oo 

Since both f{t) and g(t) = for t < 0, when we set f(t) = g(t) in the above integral equation we get, 



\G(uj)\ 2 dcu 

OO 



< I 

Also, since G(x + iy) is the Fourier Transform of exp(—yt)g(t) from Eq.(10), then by Parseval's theorem 

\G(x + iy)\ 2 dx = 2ir / e - 2yt \g(t)\ 2 dt 

-oo JO 

roc 

< 2ir / \g{t)\ 2 dt (y > 0) 

Jo 

< K (y>0) (14) 

where Eq.(13) has been used in the last step. Equation (14), expresses the property that if a causal trans- 
form is square integrable along the real axis, then once it is analytically continued into /+, it is also square 
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intcgrable along any line parallel to the real axis for < y < oo. This is a very important result and leads 
to another condition that \G(x + iy)\ — > for y > when x — ► ±oo. This also forms part of Titchmarsh's 
Theorem 93, Lemma [21], which we shall define in the upcoming section. 

Now to apply these results to our situation, first we note that we are dealing with electric fields whose intensity 
must be finite, hence both the input and output can be taken to be square integrable. Mathematically, 



f 

J — < 



\E(u))\ 2 duj < oo (15) 



which holds for both the input E^ n and output -E ut- Considering the frequency dependent terms, we note 
that, Eq.(4) is of the same form as Eq.(7). Here X{u) = E on t(S,uj), F(lu) = £7j n (0,w) and 

G(w) = e lUJn ^ 5 ' c (16) 

Hence ^ 

9 [S,t) = ^- [ e lUn{UJ)5/c e- lUT duj (17) 

where r = t — t' and 

/oo 
g(5,t-t')E in (0,t')dt' (18) 
-oo 

The integral for G(u>) is defined for negative and well as positive frequency. The refractive index is usually 
defined for positive frequency only. Since E- m and -E ut are both real, g must be real hence 

n(-uj) = n*{uj) (19) 

for real u> which extends the definition of n(u) for the negative frequencies. Due to the spatial separation 5 
between input and output we need to invoke the relativistic causality condition rather than the "primitive" 
causality condition quoted earlier. We need to restrict out output signal so that no signal can propagate 
with velocity greater than the velocity of light, c. This implies that E 0U i(S, t) can only depend on E^ n (0, t') 
when t' < t — S/c, since the output time t = t' + S/c. The signal takes a time S/c to traverse the medium. 
Hence g(S, r) = when r < 5/c. Furthermore, from Eq.(16), 

G(w) = e iUn ^ c = r g(S,T)e lU;T dT (20) 

JS/c 

Note that usually |G(w)| < 1 even without the transmission coefficient since the complex refractive index 
results in exponential decay. It is usually the case that the output signal intensity is less than or equal to 
the input signal intensity. The above equation can be written as, 



G(lu) = e W)-i]tf/c = f°° g( S , t' + 5/c)e lUJt 'dt' 
Jo 



(21) 



The right hand side of this equation has a regular analytic continuation in i+ (the upper half of the complex 
plane), this implies that the left hand side of the equation does also [20, 21]. Titchmarsh [21], gives a very 
clear criterion for analytic continuation of a function which is zero for all negative values of its argument. 
His theorem 93 has been restated by Toll [22] as; 

"A function of integrable square is zero for all negative values of its argument if and only if its 
Fourier transform is a causal transform" . 



Thus the causality condition that g(S, t' + 5/c) = for negative t' is equivalent to the requirement that G(u>) 
be a causal transform. 
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4 Titchmarsh's Theorem 



Theorems 93-95 of Titchmarsh can be restated as follows: 

Let G{oS) be an analytic and square integrable function along the real axis. (ie. Q{uS) G L 2 (— 00,00)). Then 
if G{oj) obeys one of the three conditions below, it obeys all three of them. 

(i) No output before the input, causality condition: If g(r) is the inverse Fourier transform of G(ui), then 

g(r) = (r < 0) (22) 

this means GioS) is a causal transform. 

(ii) G{x) is, for "almost all" x, the limit as y — > of an analytic function G(x+iy) that is regular (holomorphic) 
in I + (the upper half of the complex plane) and square integrable over any line parallel to the real axis. 
Mathematically this means, 

, ■ - x - 

\G{x + iy)\ 2 dy<K (y>0) (23) 



/ 



where K is a finite real constant. See our derivation of Eq.(14) above. It follows that 

lim G{x + iy) = (y>0) (24) 

x — > ±00 

The above limiting result comes from the following consideration; Cauchy's theorem for y > gives, 

If GW)dv 

yicjj = / iuj, v complex) 

2m J r v — LU 

= ±rw)*L (25 ) 

2m J ^ u-uj 

Knowing only that G {u>) is square integrable is not sufficient to imply causality. (Nussenzveig [20] , gives a 
counter example in his book page 28.) A causal transform which is square integrable along the real axis 
is also square integrable about any line parallel to the real axis in the upper half of the complex plane. 
Consider a square contour T in the use of Cauchy's theorem above, see Titchmarsh pp 125-128 [21] Theorem 
93, Lemma. (See also the treatment by Nussenzveig pages 22-24 [20]). The integrals along the vertical sides 
must tend to zero as x — > ±00. This calculation is repeated in the Appendix for completeness, 
(iii) Re G{oj) and Im G(oj) are conjugate functions and as such arc Hilbert transforms of each other, ie. 
G(u) r ) = Gr{u r ) + iGi(uJ r ) is a causal transform if and only if its real and imaginary parts satisfy; 

y r {u) r ) = — / , 

v J -00 v-u r 

(26) 

where the frequency is treated as real oj — uj r . Each of the above expressions implies the other. Here 
P indicates the principal part of the integral, (sometimes denoted by a bar through the integral sign) to 
be evaluated at v — uo r . These relations between the real and imaginary parts of G{u) are the dispersion 
relations also known as Plcmelj formulas. 



5 Dispersion Relations 

Returning to our specific example, and Eq.(21), since 5 is very small, but otherwise arbitrary, u [n(u>) — 1] 
must have the same property (as G(u>)) of analytic continuation in I + [20]. Nussenzveig also shows that 
the refractive index is regular in the upper half of the complex plane. However, this does not exclude the 
possibility of singularities on the real axis. This condition is also not sufficient to show that the real and 
imaginary parts of [n(u>) — 1] must be conjugate functions, or equivalently, that they obey dispersion relations 
and hence imply causality. We must also show that [n(uu) — 1] — > as u — > 00, which is part of condition 
(ii) in Titchmarsh's theorem above. This is in general not easy to do without resorting to a physical model 
for the refractive index. 
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5.1 Using a Model for refractive index 



The usual approach in the literature, taken by both Nussenzveig [20] and several text books [10, 11, 12], is to 
assume that for very high driving frequency, the frequency will be much higher than any binding frequency of 
the electrons in the medium, and so the electrons will behave as though they are free. The Lorentz equation 
of motion becomes 

mx = eE(uj)e- lU}t (27) 

which would involve a trial solution of the form x(t) = x(lu) exp(— iuit) and hence a result x (uo) = -e£(w)/(mw 2 ). 
Now since the polarization of the medium is given by 

V{u) = Nex{u) = X (lo)E(lu) (28) 

where TV is the number density of the electrons, m is the electron mass, and in Gaussian units, the electric 
susceptibility x(w) is, 

1 Ne 2 
X M = -[ n »-l]=— ( w -oo) (29) 

Hence, n(ui) — ► 1 as u) — > oo from physical considerations. This shows that [n 2 (o)) — l] £ L 2 {— oo, oo), is 
square integrable along the real axis. Here, x( w ) plays the role of G{lo). The causality condition, that you 
get no polarization before the electric field is applied, implies that x( w ) = for t < 0, so that condition (i) 
of Titchmarsh's theorem is satisfied. Hence, 



J — t 



\n 2 (x + iy)-l\ 2 dy < K (y > 0) (30) 



where if is a finite real constant. Thus, from physical considerations, [n(x + iy) — 1], which is analytic in 
7_l_, from the above discussion, must approach zero as x — > ±oo when y > 0. This implies, 



f 



\n(x + iy) - l\ 2 dy < K' (y > 0) (31) 



where K' is a real, finite constant. Since this inequality satisfies condition (ii), then [n(u)) — 1] must also 
satisfy the dispersion relations and be a causal transform. 



5.2 No physical model for refractive index 

We have yet to make clear that the physical model is not essential to the derivation of the dispersion relations. 
We shall now obtain the same results as in section 5.1 without invoking the physical model for refractive 
index. We begin by considering the polarization Eq. (28). This is clearly of the form Eq. (7), where the 
polarization V(oj) = X(u), E{lo) = F(uj) and x( w ) = G(u). The intensities of the electric fields must remain 
finite so it is reasonable to assume these quantities are square integrable along the real axis. x( w ) ma y also 
by square integrable or it may satisfy the weaker restriction of being bounded |x(w)| 2 < Kg, where Kq is 
a constant. In this case it is possible to form a new function by the method of subtraction (see discussions 
section 6.) which is square integrable. For physical media (of finite energy) it is reasonable to assume 
that the function x( w ) is decreasing and that it tends to zero for high frequency in all physically realizable 
situations. Electrons in a real media cannot oscillate at infinite frequency which requires infinite energy and 
thus the medium cannot become polarized at infinite frequency. Furthermore, the causality condition, that 
you get no polarization before the electric field impinges on the medium, requires that x( w ) is also causal. 
Since x( w ) is proportional to n 2 {uo) — 1 (see Eq. (29) ) then n 2 {uo) — 1 is causal, and also square integrable. 
Now it appears that n 2 (u) — 1 satisfies Titchmarsh's condition (i), so that implies it satisfies all conditions. 
We can therefore write down a dispersion relation for n 2 {ui) — 1 and in doing so we assume that (by condition 
(ii) of Tichmarsh's theorem) n 2 {uj) — 1 — > as oj — > oo. This in turn implies that n(uj) — » 1 as u> — > oo note 
without having invoked any physical model for refractive index. 

Returning to our earlier discussion, in section 3. and at the beginning of section 5., we had found that 
n(uj) — 1 could be extended into I + . We may assume that n(u) — 1 would be at least bounded. The 
remaining requirement was that n(uj) — 1 ^ as w ^ oo in order to satisfy Titchmarsh's condition (ii). In 
considering n 2 (uj) — 1 we have found that indeed n(ui) — > 1 as uj — > oo. Now we have satisfied condition (ii) 
of Titchmarsh's theorem for n(w) — 1 and thus we may write down the dispersion relations for ti(lu) — 1. 
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5.3 Calculation of the Dispersion Relations 



Now that we have established that n(uj) — 1 satisfies Titchmarsh's theorem, whether we invoke a physical 
model for refractive index or not, we may calculate the dispersion relation. Cauchy's theorem expresses 
[n(oj) — 1] at any point in I + , in terms of a contour integral, for example a large semi-circle in the upper 
half plane. Here we define frequency as a complex quantity to = u> r + iu>i- Since we have shown that the 
refractive index (for real frequency) n(u r ) — > 1 as u) r — > oo the contour integral reduces to an integral over 
the real axis, as indicated by condition (ii) in Titchmarsh's theorem above, and derived for a rectangular 
integral in the Appendix. Hence we have, 

Taking the limit as to approaches the real axis, we set complex to — > to r + it where e is a small term, the 
radius of a small semi-circle in the upper half plane about the real frequency ui r . We shall assume that ia T 
represents a pole on the real axis. We obtain, 



1 V ' J 27r»y_ 00 v-u r -% 



The denominator can be written as 



v — uj r — it \v — LJ r 
hence, in the limit where t — > we have, 



(33) 

it 



P + i7vS(u - UJ r ) (34) 



- . . - P f°° n(i/)-l ]dv . . , 

2 [n(u r ) - 1] = — / + n(w r ) - 1] 

m J„ 00 v - uj r 



J-oo (V ~ OJ r )(v + U) r ) 

Using our symmetry relation Eq.(19) we find n r (u>) is an even function and rii{io) is an odd function. We 
ignore the odd integrands which give zero upon integration. Thus we can write the real and imaginary parts 
of n(ui) where the frequency uj r = w is now purely real, as 

2P f x vim [n{y)\ dv 



ip r 

Re [n(u)} = 1 + — / 

JO 



_ 1 + £5r«62* (36) 



7T 







Im[n M ] = -_y o ^ (37) 

where a is the absorption coefficient defined earlier. These results are consistent with the findings of the 
text books [10, 11, 12] and also Nussenzveig [20], Toll [22], Landau and Lifshitz [23] and Rauch and Rohrlich 
[24]. 



6 Discussion and Conclusions 

We have shown in section 2, that although the Scharnhorst effect suggested that faster-than-c signals were 
in principle possible, in practice it would be impossible to detect any such increase in c. We further suggested 
that the original calculations of the effect may be flawed, and that the effect implies a violation of SR. Later 
sections of this paper relate to the case of physical media, by which we mean non-vacuum, finite energy 
media in non-curved spacetime. We show that physically a result n(u)) — > 1 as u — > oo is the only reliable 
result based on electron energies being less than infinity. 
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Historically, the result of Eq.(36) was first derived by Kronig (1926) [25] and the equivalent result for 
the dielectric constant was treated by Kramers in (1927) [26]. Kronig was interested in the physical model 
behind the derivation of the refractive index and how many atoms are required before you can sensibly use 
the refractive index idea for a bulk material. Kramers first employed Cauchy integrals to derive the above 
dispersion relations. However, neither paper puts any emphasis on causality. This comes in much later by 
Kronig (1942) [27]. More recent tutorial accounts of the dispersion relations which include a discussion on 
causality can be found, for example see [28, 29]. Text book accounts [10, 11, 12], will mention that n(ui) — > 1 
as lo — > oo but will only argue based on the physical model for the refractive index. We have shown above that 
if [n(u)) — 1] satisfies the dispersion relations then, regardless of how it was derived, Titchmarsh's theorem 
proves that [n(u>) — 1] obeys causality and must also therefore have the limiting result that n(u>) — ► 1 for 
lo —* oo by condition (ii) of Titchmarsh's theorem. This, to the authors knowledge, has not been stressed 
in the literature. No formal proof of n{tS) —* 1 as lo — > oo has been previously given, it has always been 
tacitly assumed from a physical model. I would like to stress that the causality condition (or the dispersion 
relations) requires this to be true regardless of the model in use. No model for refractive index can claim to 
have a different limiting value on n(u) asw^oo without violating causality for n(oo) — 1. 
Landau and Lifshitz [30], have a footnote in their book relating to the dispersion relation for the dielectric 
contant e. 

"The property e — ► 1 as lo — > oo is not important: if the limit e(oo) were other than unity, we 
should simply take (.{lo) — e(oo) in place of e(o>) — 1 ..." 

The same argument would apply to n(u>) — 1, since it is n(u>) — 1 — > as lo — > oo which is implied by 
Titchmarsh's theorem to be a condition for the dispersion relation to hold. If in a case when n(oo) ^ 1 
we used instead n(u>) — n(oo) then you would regain a valid dispersion relation in a sense that effect would 
not precede cause. It would be necessary to replace n(co) — 1 by n(co) — n(oo) everywhere in the dispersion 
relations, which is equivalent to a dispersion relation with one subtraction at n(oo), to be discussed below. 
Landau and Lifshitz were clearly not conserned with relativistic causality. The dispersion relations require 
primitive causality only, effect cannot precede cause. The author believes the above footnote to be somewhat 
misleading, especially now when physicists are conserned with signals traveling at faster-than-c velocity. The 
value n(oo) = 1 is certainly important. It is required that n(oo) = 1 to show that a front velocity of a signal 
propagates at c/1 and not c/n(oo). Pulse propagation is governed by Maxwell's equations which are certainly 
relativistically causal. Recent articles [13], on superluminal signal propagation based on unphysical models 
which allow n(w) — > (3 when lo — > oo where /? < 1 must therefore be dismissed as a violation of relativistically 
causal behavior. 

Finally, we would like to mention a paper by Weaver and Pao, [31]. They derive dispersion relations (for the 
wavenumber k(a>)) for a general class of linear homogeneous or inhomogeneous media. The proof proceeds 
without a priori knowledge of the dispersion equation for the medium, or rather without knowledge of the 
phase velocity Coo = lim (ui/k(Lu)) as lo — > oo. They represent k(u) as a Hcrglotz function. A Hcrglotz 
function k(u>) has the following properties; k(u>) is regular in I + and Im k(cu) > 0, (sec Appendix C page 
393 in Nussenzveig's book [20]). The Hcrglotz function has the desirable feature that its imaginary part 
is non-negative. The medium is assumed to be passive, no energy can be given to the wave (|G(o>)| < 1) 
which unfortunately negates the possibility of an amplifying medium, which may also be causal although its 
square integrability may be harder to ascertain. The author finds the use of Herglotz functions unnecessary 
(although mathematically valid). The required proof stems from Titchmarsh's theorem 93, Lemma, which 
Weaver and Pao seem aware of. Weaver and Pao have also gone to great lengths to explain the use of 
subtractions in the dispersion relations, also explained in detail by Nussenzveig [20]. For example, if the 
square integrability condition on G(lo) cannot be satisfied, but the weaker condition that G(lo) is bounded, 
|G(w)| 2 < Kq, can be satisfied, then we may construct a new function 



where H(lo) is square integrable and has no poles in I + and hence satisfies the dispersion relations. 



H(lo) 



G{lo) - G(lo ) 



Im loq > 



(38) 



LO — LOq 




real lo 
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where G(u>o) is called the subtraction constant. Taking the real part we find, 



Re G(oj) = Rc G(lu ) + 



/C 



Im 



G{v) - G(co ) 



V - UJ 



dv 



V — LO 



(40) 



This is known as a dispersion relation for G(ui) with one subtraction. Often, subtractions will occur at 
uj a = or ujq = oo. More than one subtraction is allowed. Using the ui = oo subtraction for n(ui) — 1 we 
find, 



Re n(ui) 



Re n(oo) + lim 

LOo— *oo 



LJ — Uq 
7T 



Im 



n(v) - njfjp) 

V - U! 



dv 



v — u> 



) — n(oo) 



v — UJ 



P f°° 

= Re n(oo) H / Im [n(v 

* J-oc 

2P f°° 

= Rc n(oo) H / [f Im n(u) - w Im n(oo)] — 



(41) 

where we have used the fact that Im n{v) is odd. This is the same as the result implied by Landau and Lif- 
shitz [23, 30], but it only requires primitive causality. For relativistic causality we also require that n(oo) = 1 
so that the front velocity of a signal travels at c/1 and not some arbitrary c/n(oo). 



We have not commented on situations which relate to highly curved spacetime. These extreme cases of 
general relativity may allow for faster-than-c signals, but we have not studied those situations here and will 
not confirm or deny any claims made in the literature. 



Appendix 



In this section wc would like to show that a causal transform of the form 



G(w) 



/•OO 

/ ff ( 
Jo 



T)e lUjT dT 



( g(j) = when r < ) 



(42) 



which is square integrable along the real axis is square integrablc along any line parallel to the real axis in 
the upper half of the complex plane, I + . This condition further implies that 



lim G(x + iy) = 

X— »±oo 



(»>0) 



(43) 



The following mathematical appendix is a summary from Titchmarsh's book [21], Theorem's 93 with Lemma 
and Theorem 95 which is related. These arguments can also be found in the book by Nussenzveig [20], page 
22. 

Consider the analytic continuation of G(u>) into I + . Let v = x + iy and let w = x + iyo and consider the 
following Cauchy integral, 

G{v)dv 



v ' 2m J r v 



(44) 

/p V — LO 

We take the contour F to be a rectangle with corners at ±i? and ±i? + iS, where both R and S shall tend 
to infinity. As we shall show, the integrals of the vertical sides of the contour are both positive and finite 
and do not cancel out. They must therefore both be zero. The magnitude of the line integral along the right 
hand side of the contour is given by, 



G(R + iy)dy 



R + iy - xq — iy 



< 



< 



max \G(R + iy)\ [ ^ — 

Jo {(R-x ) 2 + (y-yo) 2 } 1/2 



max \G(R + iy)\ In 



(S-y )+[(R-x a ) 2 + (S-y ) 2 ]- 



l(R-x )2+ : 



,1/2 



-2/o 



(45) 
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As R — > oo and S — > oo this implies that 



max + -» (46) 

0<y<S 



since the natural log term tends to ln(l + \/2) as both i? and S 1 tend to infinity. Similarly the magnitude of 
the integral on the left hand side of the contour is, 



I 



G(-R + iy)dy 



o -R + iy - x Q - iy 



G(-R + iy)dy 



o (R + x ) - i(y - y ) 



< max \G(-R + iy)\ f 

~ 0<y<S l 1 yn J Q 



dy 



< max \G(-R + iy)\\n 

~ o< y <s l y yn 



{(R + x r + (y- yo r} 1/2 
(S-y )+[(R + x ) 2 + (S- yo ) 2 ] 1/2 



[(R + x Q y + yl} 1/2 -y 



(47) 



The natural log has the same limit as above in the case when R and S* ^ oo. This leads to the same result 
as Eq.(46), now written for the limit R — > — oo. Hence, 



lim \G(R + iy)\^0 



(y > o) 



(48) 



What remains of the contour integral are the bottom and top line integrals which can be written, for R — > oo, 
as; 



According to Schwarz's inequality, 



30 G(x)dx 
oo x ~ x - iy 



G{x + iS)dx 
x + iS - x Q - iy 



(49) 



J — ( 



G(x + iS)dx 



x + iS - xq - iyo 



\G(x + iS)\ 2 dx / 

-oo J —i 



dx 



< 



-oo 

nK 



{x-x Q ) 2 + (y-y ) 2 



(50) 



where we have used Titchmarsh's theorem, condition (ii) in the last step. As S — > oo we see that the above 
expression tends to zero. Hence, 



1 f c 



G(y)dv 



v — UJ 



(Im co > 0) 



(51) 



which expresses the value of G(u) at any point in the upper half of the complex plane in terms of values on 
the real axis. 
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